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The Watson-Migdal approximation scheme to take into account final state interac- 
tions is shown to give the actual threshold momentum dependence of the reaction 
pp — ► ppn°. However, by an explicit plane wave replacement of the final state wave 
function it is stressed that not too much physical significance should be given to 
\ the proportionality coefficient extracted using this procedure. The plane wave ap- 

proximation is not physically reliable even after introducing the Watson-Migdal or 
a more sophisticated final state interaction factor, since direct production (impulse 
term) is missed. Also with short range interactions there can be discrepancies of a 
factor two. 
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During the 90 's meson production at threshold has raised much experimental 
activity and interest with cooler facilities producing data with unprecedented 
accuracy and energy resolution in a region where only a single partial wave 
amplitude should contribute. One of the exciting results is the threshold cross 
section of the reaction pp — ► ppn° [1,2], others involve e.g. 77 meson production 
[3]. This activity has also revived theoretical interest with many mechanisms 
proposed for each reaction. Discussion of these mechanisms is not the main 
purpose of this paper. 

Most theoretical approaches are based on the established and apparently in 
this case numerically well justified DWBA including various irreducible pro- 
duction mechanisms between initial and final states calculated for realistic 
interactions. There is, however, one class of works that use strong assump- 
tions about the possibility to approximate the NN wave functions either by 
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using apparently a plane wave Born approximation or including just the low 
energy final state interaction (FSI) approximately. In the latter it is assumed 
that the effect of the FSI can be factorized from the transition matrix ele- 
ments and be included after the matrix elements have been calculated using 
simple analytic wave functions, notably plane waves. Among such approaches 

- while having the value of introducing the mechanisms at a fundamentally 
more basic level as relativistically covariant or by chiral perturbation theory 

- are Refs. [4,5] on n° production and Refs. [6] on r\ production. The aim of 
this paper is to test how well such approximations reflect the numerical (and 
partly physical) reality by an explicit calculation of pp — ► ppn°. 



At low energy or threshold scattering and reactions an often used approxi- 
mation to take into account the final state interaction is the Watson-Migdal 
procedure [7]. This consists of taking first the asymptotic scattering (final) 
state wave function with the relative NN momentum pj 

Uf(r) sin(pfr + S)e lS 
p f r p s r 



(only the nucleon S-wave is considered here) and then extending this asymp- 
totic form to the range of nuclear forces (where p/r ~ 0; r cancels against the 
volume element in integration) 

u f( r ) _^ sin5e i5 = 1 
Pf pf pf cot 5 — ipf 



Here it is usual to make a further approximation in terms of the scattering 
length by replacing 

Pf cot 5 ~ , (3) 

(X 



where a is the large scattering length. So in S waves one expects the threshold 
behaviour of the cross section to be 

g ~ ^(cot^+i) xp ^rT^ xp/> (4) 



where the pf in the numerator is the momentum dependence of the phase 
space. If the infinite ranged Coulomb force is present in the final state the 
cross section should behave like 

a = -y^—, \i 7T x phase space x const. = Fpsi x phase space x const. ,(5) 
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where 5 is now the Coulomb-strong phase shift and Cq is the Coulomb pene- 
tration factor 



OJ = -5*=!-, (6) 



The remaining interaction matrix element is assumed to be relatively constant 
as compared to the fast varying Fpsi- However, it is easy to erroneously imply 
some direct specific meaning(s) for the constant extracted in this way from 
experiments. There is nothing in Ref. [7] which would relate the actual re- 
action matrix element to this constant, since the radial dependence obtained 
from the asymptotic form is meaningless at short distances. In particular this 
constant is not to be interpreted as the Born approximation to the reaction 
matrix element and this interpretation is not advocated in the original work. 
It only reflects the general final momentum or energy dependence embedded 
in Fpsi for small momenta pf. 



To exemplify this in a specific reaction with the strongly attractive nearly 
bound 1 S'o two-nucleon final state I discuss very explicitly the reaction pp — > 
ppn close to threshold. This has both complications: the strong interaction 
and the Coulomb interaction in the final state. The energy discussed is 290.7 
MeV corresponding to pf (max) = 0.35 fm _1 and is about 5 MeV above thresh- 
old in the c.m.s.. The formalism has been given elsewhere [8] and will not be 
repeated here as not very essential for the present argument. Suffice it to say 
that the pion is produced primarily by the Galilean invariant ttNN vertex for 
each nucleon obtained from the pseudovector pion coupling 

/ - f=.=. T 



H^NN = — O- \ q T ■ 0- — [p T ■ <j) + T ■ (f) p \ . (7) 



Also it is necessary to include at least pion s-wave rescattering from the second 
nucleon and some additional mechanism, such as a heavy meson exchange with 
a nucleon z-diagram in the production vertex to fit the data (without this, 
direct production is a gross underestimate as compared to data). In reactions 
where the NN isospin changes this s-wave rescattering even dominates at 
threshold, but in the present case it is smaller than the direct axial charge 
term. To avoid unnecessary complications the treatment here considers only 
nucleons, no A components. 

The conventional mechanisms lead to the exact transition amplitude at the 
two-nucleon level 



< 1 S |tf pro d| 3 Po > = ^=-L- {(1 - q I «5Kr)ji(9r/2)«i(r)dr 

PfPi v /uJ^m 7T + I 2M J J 
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+ ^Ju' f *(r)U q r/2) Ul (r)dr 

(2 - ^) / u}(r)f(r)j ( q ^) Ui (r)dr 



^ / <(r)/(r)io(f ) Ul (r)dr 



}• (8) 



Here the first term arises from the standard p-wave nN coupling proportional 
to the pion momentum a ■ q and should be small at threshold due to the 
presence of the pion wave function ji(qr/2) in the integrand and the factor q. 
The Galilean invariance term does not have this suppression, just the factor 
ujqjM < 1. The last two terms involve s-wave pion rescattering off the second 
nucleon with the isospin symmetric amplitude Ai. The Yukawa function f(r) 
and its derivative f'(r) arise from the propagator of the intermediate pion and 
have a range rather similar to OPE [9]. The wave function derivative in the 
present normalization (see Eq. (1)) is -u'(r) = r d(u(r) / r) / dr . One may note 
that the latter terms involve an "interaction" of finite range, while the first 
two overlaps have an infinite range. It is therefore of particular interest to see 
how well short range expansion of the final state wave function, the basis of 
the Watson-Migdal method, fares here as compared with an exact calculation. 
Of course, correlations are necessary for the reaction to happen at all, since an 
on-shell nucleon cannot emit a real pion. So, actually the reaction does have 
a finite range, though the integrals to be evaluated extend to infinity and no 
convergence factor is present. This expectation will be borne out in numerical 
results, where little difference is found between the energy dependences of the 
first and last terms. 

The above integrals can be calculated numerically [8]. The dominant terms 
near threshold are the Galilean invariance term (the second) and the rescat- 
tering terms. When the pion has enough energy, the first term becomes com- 
parable (in fact, physically it will need also the inclusion of the A). For the 
present study of the effect of the final state interaction, the dependence on the 
relative energy (or momentum pj) of the final state nucleons is particularly 
relevant. So special care with also the long ranged Coulomb force is impor- 
tant, since it has large influence at low energies. In all cases presented here 
the initial state is the same distorted 3 P wave function calculated from the 
Reid soft core potential [10]. This potential is also used for the final state in 
the "exact" case. 

Fig. 1 shows the ratio of the squared transition matrix element (calculated 
using the exact wave function) to the final state interaction factor F FS i as a 
function of the final state relative pp momentum with the inclusion of differ- 
ent mechanisms. The dotted curve shows the small effect of the direct axial 
current coupling (nongalilean first term, multiplied by 10) f/m 7T + a ■ q, where 
q is the pion momentum. The dashed curve includes also the much more im- 
portant axial charge coupling (Galilean part) —f/m n +u q (r- (p + ff)/2M. In 
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Fig. 1. The ratio of the squared matrix elements and the final state interaction factor 
Fpsi defined in Eq. (5). Dotted curve: direct axial current contribution (nongalilean) 
multiplied by 10; Dashed: full direct Galilean production operator H^nn used; Solid: 
full result including also pion s-wave rescattering; Dash-dot: full result with added 
HME to fit the data at threshold (divided by 10). 

the solid curve also pion s-wave rescattering is added to these impulse terms 
- in several pionic reactions this is sufficient and satisfactory. By far, most of 
this contributes through the nongalilean part. The minute direct nongalilean 
result has a qualitatively different energy dependence, since it has to vanish, 
when q -> 0, i.e. when p f -> p fymax . 

However, this reaction needs some additional mechanism [1] and in the dash- 
dot curve a phenomenological heavy meson exchange is added to reproduce 
the experimental cross section [11]. It can be seen that in this momentum 
range, apart from the FSI factor, the momentum dependence is moderate for 
all mechanisms and the validity of the Watson-Migdal procedure is confirmed. 
It may also be noted that even the "long-ranged" galilean impulse term does 
have the same momentum dependence as those with an explicit short-range 
interaction. This is, of course, due to the fact that the correlations necessary 
for the reaction are generated by meson exchanges. If the Coulomb interaction 
is removed, the results are nearly the same: the ratios increase by 10-20% and 
the energy dependence is marginally weaker. 

As another step we try to interpret the matrix element as arising simply 
from a use of a plane wave for the final state, trying to take the FSI into 
account just by the factor Fpsi- Fig. 2a shows the behaviour of the squared 
transition matrix for the direct (impulse) contribution as a function of the 
relative final two-nucleon momentum pf. The solid curve employs the exact 
final state wave function (likewise in the other figures). To make a "Born 
approximation" with respect to the final state one replaces Uf(r)/(j)fr) in Eq. 
(8) by joipjr) = sin(pjr)/(pjr). The resulting negligibly small approximate 
results are not shown in the figure. 

In fact, both approximate impulse amplitudes are much smaller than the cor- 
responding exact ones, the nongalilean one by two orders of magnitude, the 
galilean one by three orders. This smallness is apparently due to the overall 
orthogonality of the three relevant Bessel functions [12] and to the fact that 
the 3 Po initial wave function does not deviate much from the plane wave (at 
290.7 MeV the calculated phase shift is —9.5°). It may be noted that also the 
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Fig. 2. Momentum dependence of the squared transition matrix. Curves: Solid) 
exact final state; Dashed) final state approximated by a plane wave and multiplied 
by Fpsi of Eq. (5) (fm 9 ); Dotted) exact s-wave rescattering without the impulse 
term; Dash-dot) the FSI factor from Eq. (9) used. Windows: a) direct ("impulse") 
terms only, b) pion s-wave rescattering included, c) also phenomenological heavy 
meson exchange added, d) only HME. In a (direct terms) the approximate results 
are negligible (smaller by 4-6 orders of magnitude). 

derivative of the final S'-wave as denned above is a Bessel function. Making 
the NN potential stronger increased these amplitudes; making it zero further 
reduced the direct amplitudes by another two orders of magnitude acting as 
a check on numerics of these infinite-ranged oscillatory integrals. 

The strong S-w&ve final state interaction removes this orthogonality making 
the amplitude sizable. These correlations influence particularly much on the 
derivative, making the matrix element of the galilean term physically impor- 
tant. Also inclusion of an explicit interaction such as pion rescattering removes 
the orthogonality in the integrand. Then apparently at threshold the approx- 
imate integral is relatively insensitive to pf, since joG°/ r ) ~ 1 f° r small r. 

With a finite ranged interaction the energy dependences of the approximate 
and exact results are qualitatively similar (dashed curves vs. solid in Figs. 2b- 
d), once the FSI factor is applied, as would be obvious from the constancy of 
the approximate integral. However, physically in the exact result it is necessary 
to take into account the significant galilean impulse term. In the galilean case 
(Fig. 2b) the impulse term is constructive with rescattering and the exact cross 
section is larger than the approximate FSI result, while in the nongalilean 
result the situation would be reversed. Even the FSI factor cannot reproduce 
the full height of the exact result, even though the approximate rescattering is 
more than twice as large as the exact (dotted curve in 2a). This shows clearly 
the physical importance of the galilean impulse term. 

Finally in Fig. 2c the overestimation of short ranged mechanisms by the ap- 
proximate factor (shown explicitly in 2d) reproduces the exact result (with a 
slight overestimate). However, this good agreement is achieved by overestima- 
tion of the short-ranged interactions and omission of direct production. 

Now it is important to note that, since also the approximate final wave func- 
tion includes the 1/pf, it is not consistent to apply the FSI factor denned 
by Eqs. (2-5) where the same 1/pf appears - the result would not be even 
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dimensionally correct. If this double counting is avoided, the resulting energy 
dependence clearly would disagree with the exact result - the dashed curves 
should be multiplied by pj. This exercise is done, because a more sophisticated 
treatment is not numerically very different (if pf is given in fm^ 1 ) and because 
it is possible that the simpler method may have been applied some times in 
the past. Also, it was seen that the energy dependence comes out correctly. 

An enhancement factor of the amplitude having a correct high-energy limit 
is based on properties of the Jost function and presented in Ref. [13] for a 
short-ranged interaction applying the effective range expansion to order p 2 

1 (pj + a 2 )r /2 



f(-pf) l/a + r pj/2-ip f 
with 



a 



= (l + yjl + 2r /a)/r . (10) 



Except for the numerator, to this order this is equivalent to the simpler form 
and is dimensionless. Also, if the unit of length is the femtometer, the numer- 
ator is indeed of the order of unity for small pf. In Ref. [4] this is used with 
a Coulomb modification on the scattering length a = 7.8243 fm and effective 
range r = 2.7058 fm given now as [14] 

l/a c =[l/a-2p fV h(r ] )]/C 2 (11) 
r 0c = r /C 2 (12) 

with the function h(r)) given on p. 263 of Ref. [13]. However, since a c and r 0c 
are not constant with pf, it may not be clear that the assumptions on which 
Eq. (9) is based are valid in the presence of the long-ranged Coulomb force. 
The numerical error arising from this complication is likely quite small but 
not under control. 

Now using this improved final state interaction factor the dash-dot curves in 
Fig. 2 show similar good qualitative energy dependence with the exact result as 
the dashed ones - only slightly weaker. However, although the normalization 
is much improved (as can be seen from the dash-dot vs. dotted curve in 2b 
and dash-dot vs. solid in 2d), still a mere multiplicative factor cannot correct 
for the missing direct production strength (the dash-dot vs solid curves in 2b 
and 2c). Consequently, even after the inclusion of the heavy meson exchange 
the total results remain now as underestimates shown in Fig. 2c, while only 
this short range contribution is presented in Fig. 2d to see the effect of the 
range alone. In this case the approximation is a slightly larger overestimate, 
about a factor of three for the simple factor. 
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Fig. 3. As Fig. 2 but without the final state Coulomb force 



Above, the Coulomb interaction has been included in the final state interac- 
tion. It is also interesting to see how the approximation fares without this 
complication. After all, Eqs. (9-10) are deduced for the uncharged effective 
range parametrization. One might consider this to describe e.g. nn — > nnn° 
(hardly physically measurable), np — > nnn + or np — > npn reactions. How- 
ever, to avoid superficial differences the kinematics is kept the same relevant to 
pp — > ppn°. Even so, the scattering length and effective range are taken to be 
the same experimental np singlet parameters a np = 23.715 fm and r 0np = 2.73 
fm as in Ref. [4]. Use of the Coulomb modified quantities above relevant to pp 
scattering would hardly be physically meaningful and would give an incorrect 
view of the effect. 

The results are shown in Fig. 3 in the same way as in Fig. 2. It can be seen 
that now both approximations are close to each other, since the numerator 
of Eq. (9) is (incidentally) close to unity, if the momentum is given in fm -1 . 
Now both approximate methods give an overestimate by factor two for matrix 
elements involving potential ranges. The total result now agrees incidentally 
with the exact result but due to the overestimate of HME and omission of 
direct production as earlier for the simpler procedure. Also, it can be seen 
that for small values of pf the functional dependence changes qualitatively. 

In summary, it was seen that the Watson-Migdal conjecture that at low mo- 
menta the reaction matrix momentum dependence can be obtained from the 
asymptotic scattering wave function (essentially the phase shift) is true in this 
reaction. However, using this procedure in the other direction by calculating 
the reaction matrix by plane wave functions and simply applying the final 
state interaction factor of the Watson-Migdal method or by the method given 
in Ref. [13] is risky and may lead to physically incorrect conclusions about 
the reaction mechanisms. At its simplest, the energy dependent enhancement 
factor Fpsi cannot be correctly used with the Born approximation potential 
matrix elements, whereas the more sophisticated procedure gives them reason- 
ably well. Both approximations overestimate short-ranged mechanisms when 
the Coulomb force is not present in the final state. 

In particular, in pp — > ppir the multiplicative procedures miss the important 
direct production mechanism. In principle, it is possible to incorporate most 
important meson exchanges after the pion production vertex. However, this 
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may lead to strong violation of unitarity and also obscures the role of the 
FSI. Inclusion of any FSI factor would then risk doubly counting parts of this 
interaction. 
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